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1. Introduction 

Performing multiloop calculations is usually most convenient within the framework of 
dimensional regularization. However, when chiral fermions are considered, one often en- 
counters technical difficulties with extending 75 to D = 4 — 2e dimensions. These difficulties 
are related to Dirac traces like 

T^,^pa = Tr(7„7^7,,7p7^7"75). (1) 

Such a trace is not uniquely defined in D dimensions, so long as one assumes that 75 anti- 
commutes with all the other gamma matrices. One can easily check that the totally anti- 
symmetric part of T^upa equals to 

T[^upa] = -D Tr(7[^7^7p7^]75) (2) 

when cychcity of the trace and {70,75} = is used, or to 

T[^^,ypa] = {D-8) Tr(7[^7^7p7<,]75) (3) 

when the two contracted 7a's are brought to each other through the whole string of other 
gamma matrices with use of {7^,7/3} = 2gai3. 

These two results agree with each other only when D = A or when the whole trace 
vanishes. Neither ambiguous results for D A nor the vanishing trace can be accepted in 
the Dirac algebra in a consistent dimensional regularization scheme. The Dirac algebra in 
an acceptable scheme must satisfy the following requirements: 

• (i) It should give unique results independently of whether some diagram is considered 
as a subdiagram, and independently of the order in which subdiagrams are calculated. 
Otherwise sub divergences could not be properly subtracted in multiloop diagrams. 

• (ii) When D ^ A, its rules must analytically tend to the rules of 4-dimensional Dirac 
algebra. 

Several schemes satisfying these requirements have been proposed What they 

have in common is that their use requires complicated or tedious algebraic manipulations. 

Calculations become much simpler if occurrence of traces with 75 can be avoided in a 
given problem. Then, one can define the Z)-dimensional 75 as 

(D-l)(D-2) „ , 



With such an exphcit definition at hand, the requirement (i) is automatically satisfied. The 
requirement (ii) would not be satisfied for diagrams containing traces with 75. (The trace 
(1^) would then identically vanish for D ^ A). However, in the absence of traces with 75, the 
requirement (ii) is satisfied, too. The matrix defined in eqn. (^) anticommutes with all the 
remaining gamma matrices, and its square is equal to unity. 

Flavor changing |AF| = 1 {F = S,C or B) processes in the Standard Model (SM) 
are described by Feynman diagrams with no traces containing 75. This statement is true 
only at the leading order in weak interactions, but to all orders in QCD and QED. While 
higher orders in weak interactions are usually negligible, one is often interested in resumming 
large logarithms ln(M^/m^^^^J from all orders of the QCD perturbation series. This is most 
conveniently done with use of an effective hamiltonian built out of light fields only. It is 
possible to define this hamiltonian in such a way that no traces with 75 are introduced. 

In the standard approach to |AF| = 1 processes |^, the applied form of the effective 
hamiltonian induced appearance of traces with 75. It was harmless at one loop, but caused 
many difficulties in the next-to-leading order (NLO) 2- loop calculations p|, |lT], |12|, \T^ . 
Applying the same scheme in even higher order computations would be very inefficient. 
Instead, it is better to start from the outset with a modified form of the effective hamiltonian 
which is free of the 75 problem.^ 

In the present paper, we present our calculation of the next-to-leading QCD corrections 
to the |AF| = 1 nonleptonic effective hamiltonian. We use an operator basis in which 
traces with 75 do not occur. Consequently, we are allowed to use fully anticommuting 75 
to all orders in QCD. Such an approach is much more natural and useful than using the 
"standard" operator basis [P, ^ In our case, the two-loop calculation can be made 
completely automatic and almost trivial. We do not face subtleties related to the use of 
Fierz symmetry arguments in D dimensions, as it was the case when fully anticommuting 
75 was used in the "standard" operator basis. The calculation is much simpler than in 
the 't Hooft-Veltman scheme [Q for 75, too. In this latter scheme, dimensional Lorentz 
invariance is violated, which causes appearance of a large number of complicated evanescent 
operators. 

Our main motivation for undertaking the present work was that its outcome was directly 



applicable in our 3-loop calculation of NLO corrections to b ^ Xs7 decay uB, 16 1. However, 



^ A similar issue concerning renormalization of four-quark operators has been discussed in ref. 
the context of QCD sum rules. 
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we find it useful to discuss tlie nonleptonic liamiltonian separately from 3-loop b Xs'y 
teclinicalities. 



Our calculation proceeds along the standard lines |T0|, 0. We evaluate tree- level 
and one-loop matching conditions for the operators in the effective hamiltonian, as well 
as their one- and two-loop anomalous dimensions. We only use another basis of effective 
operators. Even with our results at hand, it is quite nontrivial to verify whether they agree 
with the previous calculations |10|, |rT[ Such a comparison is essential because of the 
phenomenological relevance of the considered quantities. Therefore, we devote a sizable part 
of the present article to show how this verification is performed. 

Our paper is organized as follows. In the next section, the effective hamiltonian is in- 
troduced. In section 3, we calculate its Wilson coefficients and the necessary anomalous 
dimension matrix. Section 4 is devoted to showing how our results can be transformed 
to the previously used renormalization scheme, and to verifying that we can indeed con- 
firm results of the previous calculations. The appendix contains a list of the nonphysical 
counterterms relevant in our calculation. 



2. The effective hamiltonian 

For definiteness, we shall consider here the AB = —AS = 1 nonleptonic effective hamil- 
tonianj^ neglecting the small CKM matrix element Vub- Generalization to other |AF| = 1 
processes and/or to processes where Vub cannot be neglected is straightforward. 

The effective hamiltonian is built out of operators Pi multiplied by their Wilson coeffi- 
cients Ci 

'^eff = T^ytlVtb'^CilPi + (counterterms) i] , (5) 

where Gp is the Fermi coupling. 

The specific structure of the operators Pi is determined from the requirement that the 

hamiltonian reproduces the Standard Model AB = —AS = 1 nonleptonic amplitudes at 

the leading order in (external momenta) /M^i/ and in the electroweak gauge couplings, but 

to all orders in strong interactions. If this requirement was applied to off-shell amplitudes, 

then EOM-vanishing operators, i.e. operators which vanish by the QCD equations of motion 

would need to be present among Pj. Here, we do not include such operators. Consequently, 

our Heff reproduces only on-shell SM amplitudes (both partonic and hadronic ones) |[T8[] . 
The hamiltonian is assumed to conserve flavors other than B and S. 
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However, we still require renormalizability of IPI off-shell Green's functions generated by 
Tie//- Therefore, some EOM- vanishing operators are present among counterterms. 

The above requirements imply that the set of operators {Pi} must be closed under QCD 
renormalization, up to counterterms proportional to nonphysical operators, i.e. operators 
whose renormalized matrix elements between physical states vanish. Such nonphysical op- 
erators can be either EOM-vanishing or evanescent (i.e. algebraically vanishing in four 
dimensions) [0, Potentially, one could also encounter BRS-exact operators (i.e. BRS- 



variations of some other operators) as nonphysical counterterms ||T9|, However, they turn 
out to be unnecessary up to three loops for the operators Pi considered below. 

A set of operators Pi which satisfies the imposed requirements consists of dimension-six 
operators Pi, Pg 

P2= {sLlt^CL){cLl''bL), 

ib) 

P,= {sLl^ir%L)Y.,{qYT-q). 

P5 = (sL7m7M2 7M3&L)E,(g7^^7^^7'^^g), 

and one dimension-five operator proportional to^ 

{si^a^'^T%R)G%, (7) 

called the chromomagnetic moment operator. Here, is the SU{?>)coior generator, and G^^, 
stands for the gluon stress-energy tensor. 

Certainly, the above choice of the operator basis is not unique. One could choose other 
linear combinations of the operators and/or redefine them by adding some nonphysical ones. 

Let us briefiy argue why all the operators given in eqns. (P) and (0) must be present in 
T^e//- The operator P2 must be present because it enters the matching condition shown in 
fig. |1]. In this figure, the VT-boson exchange is approximated by the effective Fermi weak 
interaction of quarks. 

All the other operators in eqns. (|^) and (|^) are generated from P2 in the process of QCD 

renormalization of various amplitudes: The diagrams with P2- vertices shown in fig. ^ require 

counterterms proportional to Pi, P2 and to the nonphysical evanescent operator e'^\ (The 

relevant nonphysical operators are collected together in the appendix.) All the remaining 
* The s-quark mass is neglected here. Therefore, sr decouples from flavor-changing processes. 




Figure 1: Tree-level matching condition. 




Figure 2: One-loop current-current diagrams. 



divergent 1-loop IPI diagrams with P2-insertion are shown in fig. ^ They are renormalized 
by a single counterterm proportional to 

(^lT%6z.)I^.G'"'^'' (8) 

where Di, is the QCD covariant derivative. For convenience in the future two-loop calcu- 
lation, we write the above operator as a linear combination of P4 and the EOM-vanishing 
nonphysical operator A^^i (see the appendix) 

iVi = -(sLT''7^6L)/^.G"^^^ + P4. (9) 




Figure 3: One-loop penguin diagrams with no Dirac trace. 



Inserting P4 into the diagrams of fig. ^, one finds that counterterms proportional to P5 
and Pg are necessary. Finally, the same diagrams with insertions of Pg require a counterterm 
proportional to P3. 

The set of operators {Pi,...,P6} is closed under one-loop QCD renormalization, up to 
nonphysical counterterms. A set which closes to all orders in QCD is obtained by including 
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only one more operator, i.e. the chromomagnetic moment operator (|^). The fact that no 
more physical operators of dimension 5 and 6 need to be present in T^e// can be shown 



analogously to ref. p2| by writing all the possible AB = —AS = 1 operators allowed by 
gauge symmetry and reducing them by the equations of motion. 

In the standard approach P, |10|, , triple products of gamma matrices occuring in P5 
and Pg have been reduced with help of the four-dimensional identity 

1f^7u7p 9,^ulp + Qvpli, - Qpplu + if^i.upal'^lb- (10) 

Doing this step in a two-loop calculation requires introducing several more evanescent oper- 
ators and leads to problematic traces with 75. Therefore, we do not perform this step, and 
just leave the operators P5 and Pe as they stand.0 We also find it convenient to leave their 
color structure as it is, rather than to remove the generators T"" with help of the identity 

In the remaining part of this paper, we shall concentrate on finding the Wilson coef- 
ficients and the anomalous dimension matrix of the operators Pi,...,P6, leaving aside the 
chromomagnetic moment operator (^). The latter dimension-five operator does not require 
counterterms proportional to dimension-six operators. Thus, it does not affect renormaliza- 
tion group evolution of their Wilson coefficients. Consequently, the operators Pi, ...,Pq can 
be considered separately, which makes the necessary formulae more compact. We discuss 
the chromomagnetic moment operator together with its photonic counterpart up to three 
loops in other publications |jT5|, |16 . 



3. Matching conditions for the Wilson coefficients and the anomalous dimension 
matrix 

We now turn to the values of the coefficients Ci, Cg. They are found by matching the 
SM and effective theory amplitudes perturbatively in a^. When performing the matching, 
one sets the renormalization scale fi close to M]y, in order to reduce size of ln(Miy//i) which 
could otherwise worsen the perturbative expansion. In the following, we use the MS scheme 
with fully anticommuting 75, and take /i = M]^ as the matching scale. 

The tree-level matching condition is shown in fig. |I]. The one loop matching conditions 

for the coefficients Ci and C2 are found by requiring that the one-loop SM diagrams in fig. ^ 

^ An essentially identical procedure was used in the calculation of one-loop corrections to the coefficient 
functions of four-quark operators appearing in the operator product expansion of two quark currents . 
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be reproduced by the effective theory diagrams in fig. The one-loop contribution to the 
coefficient C4 can be found by on-shell matching of the IPR amphtudes presented in fig. |. 
The coefficients C3, C5 and acquire no one-loop contributions. 
The obtained coefficients Cj(/i = M^r) read 



Ci{Mw) 
C2{Mw) 



+ 0{a 



A-n 

1 + 0{al 



4:71 



E{x) - - 



+ 0{ai 



(11) 
(12) 

(13) 

(14) 



where 



E(x) 



x(18- lla;-a;2) ^^(IS - 16a; + 4^^) , 2, 

H ^-77z ^lna;--lna;. 



12(1 -a;)3 



6(1 



and X = m^/M^. 



Figure 4: One-loop SM current-current diagrams. 




(15) 





Figure 5: One-loop matching for C4. Dots stand for 0{1/M^), a counterterm, and for a 
diagram with insertion of the chromomagnetic moment operator given in eqn. (0). 



Many flavor changing processes occur at energy scales mught << My/. Perturbative 
expansion of their amplitudes is efficient when the renormalization scale yU is close to mught 
rather than to My/. Wilson coefficients Cj at /i ~ mught are found from Ci{Mw) with help 
of the Renormalization Group Equations (RGE) 



d 



(16) 



^ Tree diagrams with counterterms need to be added on both the SM and the effective theory sides before 
the matching is performed. 
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Here, 7(/i) is the anomalous dimension matrix which has the following perturbative expansion 

m = + #^^7^^^ + 0{al). (17) 



The anomalous dimension matrix in the MS or MS schemes is found from one- and 
two-loop counterterms in the effective theory, according to the following relations: 

7(0) = 2a", (18) 
7(1) = Ad^^-2bc. (19) 

The matrices a^^, a^^ and b in the above equations parameterize the MS'-scheme counterterms 
proportional to Ci, ...,Cq in eqn. (||) 



(counterterms)i 



.k=l k=l 



O^g / I 22 12 \ o 



(4vr) 



k=l 



k 



+ (two-loop chromomagnetic counterterm) 

-|- (two-loop nonphysical counterterms) -|- 0{a^). (20) 

The matrix c occurs in the one-loop matrix elements of evanescent operators. Let us 
denote by {E\}^)iioop any one-loop amplitude with an insertion of some evanescent operator 
Dependently of what external lines are chosen, such an amplitude can be given by 
the diagrams in fig. ^ or in figs. |^ and ^. Pole parts of such amplitudes are proportional to 
evanescent operators. The remaining parts in the limit D — > 4 equal to linear combinations 
of tree-level matrix elements of the physical operators and A^i. 
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{Ei^)lloop = — ~ \dkj{El ^)tree + ekj{E\ '^)tree ~ Cki{Pi)tree 

j=0 ^ 1=0 

+ Xk{{sLa^'''T%R)G%)tree ' nk{N^)tree + 0{e). (21) 



The coefficients in these linear combinations (which define the matrix c) are independent of 
what particular external lines are chosen. However, recovering the matrix c usually requires 
considering several choices of the external lines, so that all the potentially relevant operators 
can enter into eqn. (^Tj) with nonvanishing tree-level matrix elements. 



In ref. pO[, one can find a more detailed explanation why evanescent operators affect the 
anomalous dimension matrix precisely in the manner we have just described. 

The matrices a^^, a^^, h and c which determine 7'-°^ and 7^^-' are found by evaluating 
various one- and two- loop Feynman diagrams with insertions of the operators Pi and E^\ 
The one-loop diagrams are shown in figs. |^, ^ and All the two-loop diagrams we calculate 



Figure 6: One- loop penguin diagrams with Dirac traces. The traces contain no 75. 



(where only Pi need to be inserted) are identical to IPI parts of the two-loop diagrams shown 
in the figures of ref. |T^. Most of them can be obtained from the one-loop ones by adding 
an additional internal gluon. 

We evaluate pole parts of these diagrams using the method we have described in refs. [0, 
^ where a common mass parameter was used as an infrared regulator. Next, we decompose 
the obtained ultraviolet divergences into terms proportional to the physical and nonphys- 
ical operators. This gives us the necessary renormalization constants and, in the end, the 
anomalous dimension matrix for the operators Pi, Pg- 



We find (both in the MS and MS schemes) 
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and 
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77600 1264 r 
243 81 J 


28808 1 164 r 
243 ~^ 81 J 


20324 1 400 r 
243 ~^ 81 J 


21211 1 622 r 
162 ~^ 27 -I 



(23) 

where / is the number of active fiavors (equal to 5 for fi G [rrib, Mw])- 

In the end of this section, we shall use the obtained anomalous dimension matrix to find 
explicit expressions for the Wilson coefficients 



a(/i) = c^r(/i) + 



(0)... , ^C«(/.) + 0(a^) 
9 



(24) 



at the renormalization scale ~ mi, which is appropriate for studying b-quark decays. 
Using the general solution of the next-to-leading RGE (given eg. in ref. we find 



where t] = as{Mw)/c(s{f^b) and 
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B' 
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1.9737 

-0.5409 
2.2203 
0.0400 

-0.2614 
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0.0971 
-0.0162 
-0.0243 
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0.0237 
0.0094 
0.0108 
1.0951 





0.0595 
-0.0173 
-0.0100 

0.0163 







-0.0218 
-0.1336 
0.0010 
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0.0335 
-0.0316 
-0.0017 

0.0023 
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(25) 
(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



For as{Mz) = 0.118 and /x = 5 GeV, the ratio r] at NLO is r] 
0.1203/0.2117- 0.568, and 



as{Mw)/as{Hb) 



-0.480 1.023 -0.0045 
12.12 -0.965 0.0230 
-0.276 1.007 -0.0041 



-0.0640 0.00043 0.00091 
-0.779 -0.0093 -0.0083 
-0.0771 0.00028 0.00077 



(32) 



where mt = 175 GeV has been used. 



10 



When processes taking place at energy scales much lower than rrib are considered, the 
6-quark must be decoupled at /i ~ rrib. The same refers to the c-quark at /i ~ rric in the 
case of kaon mixing or decays. The number of active flavors / in the anomalous dimension 
matrix changes at the heavy quark thresholds. One-loop matching between effec- 

tive theories needs to be performed then (see eg. ref. p6|). Explicit formulae for the Wilson 
coefficients become more complex in such cases. 

Having evolved the effective hamiltonian to a low-energy scale, one calculates its matrix 
element between some physical states of interest. The actual way of calculating the matrix 
element depends on the process under consideration. One must always make sure that the 
same renormalization scheme is used in the renormalization group evolution and in evaluating 
the matrix element. 

We shall not consider any examples of calculating matrix elements here. Instead, in the 
next section, we shall demonstrate that our anomalous dimension matrix is in agreement with 
the one found previously in the standard operator basis. Consequently, phenomenological 
results obtained in our basis cannot be different from what is already known. The advantage 
of introducing our new basis is not improving the existing two-loop phenomenology, but 
rather forming a much more convenient starting point for even higher loop computations 



like the one in refs. [15, Id 



4. Transformation of the anomalous dimension matrix to the "standard" basis 

The present section is devoted to demonstrating that the anomalous dimension matrix 
(|22|) -(p3|) found by us in the new operator basis © agrees with the one found previously in 
the following "standard" basis of operators ||, |10|, plj] : 

02= {sll.cDic^rbl), 
03= {sll,hl)Y.,{qh<). 

06= {siiA)^MriR)^ 
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supplemented with the following evanescent operators: 

Of = {ni.^i.^i^AMr'r'r'bi) + (-le + 4e) o^, 

Of = (s27m7M.7M3c2)(ci7^^7'^^7^^&i) + (-16 + 4e) O2, 

Of = (^27m7M.7M3&2) E,(gf 7'^^7^^7'^^gf ) + (-16 + 46) O3, ^^^^ 

Of = (s27m7M.7M3&^) E,(gf 7'^^7^^7'^^gg) + (-16 + 4e) O4, 



Of = (s27«7m.7m3&l) E,(g^7^^7^^7'^^g^) + (-4 - 4e) O5, 
Of = (s27«7m.7m3&l) E,(g^7^^7^^7'^^gR) + (-4 - 4e) Og. 



:.( 

Such evanescent operators are necessary to renormalize one-loop current-current diagrams 
(fig. ^ with insertions of the operators Oi,...,Oq. We have found their explicit form by 
imposing the requirements given in eqns. (4.3) and (4.8) of ref. JlOl up to O(e^). 

We shall pass from the new operator basis (p|j6^) to the old one ( |33| , p^ in a series of 
subsequent redefinitions of operators. First, we shall redefine the physical operators by 
adding some evanescent ones to them. Next, some of the evanescent operators will be 
redefined by adding "e x (physical operators)" to them. The final (and easiest) step will be 
an e-independent linear transformation of operators. 

Before going into details, let us mention how the renormalization constants and anomalous 
dimensions transform when passing from one operator basis to another. Let us start with a 
set of physical operators {Pi} and a set of evanescent ones {Ej^^}. The quantities we would 
like to consider are the one- and two-loop anomalous dimensions 7*^'^-* and 7^^-', the one- loop 



renormalization constants b and d, as well as the matrix c (see eqns. (|y) and (21)). 

Suppose we want to pass to a new (primed) operator basis which differs from the initial 
one by adding some linear combinations of evanescent operators to the physical ones 

p:=p^+J:w.kEl}\ Ei'^ = E^^\ (35) 

k 

The quantities under consideration transform as follows: 

= ^(0), ^'(1) = ^(1) - \Wc, 7^°)] - 2/5o#c, (36) 

6' = 6 + l^Wiy, c' = c, d' = d. (37) 

where /^o = H — §/ is the well-known one-loop /3-function coefficient in QCD. 

Next, we redefine the evanescent operators by adding "e x (physical operators)" to them 

P: = Pi K^'^ = E't'^ + eY^Uk^Pi (38) 
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The corresponding transformations read 

y'(o) = ^/(o)^ ^ ^/(i) ^ pf>^ ^/(o)j ^ 2/?o6'?7, (39) 

b" = b\ c" = c' + ^U^'^'^ ~d'U, d" = d'. (40) 

Finally, we perform e-independent linear transformations of the physical and evanescent 
operators (separately), passing to the triply primed basis 

P'' = E Ef'^ = E M,iE';^'l (41) 

j I 

The transformations are now rather trivial 
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'"W - i?7"(°)i?-\ 7"'(^) = ^7"(i)i?-\ (42) 



h"' = Rh"M-\ c"' = Mc"R-\ d"' = Md"M-\ (43) 

Before we start to transform our particular basis of the physical operators Pi (eqn. (|^)), 
and the one- loop evanescent ones (eqn. (§^)), we extend it by introducing four extra evanes- 
cent operators 

E'i^ = P5-10P3 + 6P3, 

E'i^ = Pe- IOP4 + 6P4, 

41)= P5_10P3 + 6P3, ^^^^ 

Ei^^ = Pe - IOP4 + 6P4, 

where tilde over an operator denotes inserting 75 in front of the last quark field. If some oper- 
ator Pi is a linear combination of terms (sLr6/,)(gr'g), then Pj is the same linear combination 
of (sir6j^)(gr'75g), where F and F' stand for arbitrary Dirac matrices. 

The operators E^^\ Eg^^ are not needed as counterterms in the initial basis. However, 
some linear combinations of them will become parts of either the physical or the evanescent 
operators in the triply primed basis which will be equivalent to the "standard" one (eqns. (|33|) 
and (§5)). 

The anomalous dimension matrices 7^°-' and 7*^^^ in the extended basis are just the same 
as in eqns. ( ^2]) and ([23|). The matrices b, c and d in the extended basis are found from 
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one- loop diagrams (figs. ^ ^ and ||) with insertions of Pi, ...,Pq and E 



(1) 
1 ) 



Ei^\ We find 
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(45) 
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? 

? 








8960 
3 

4480 
9 

320 
3 

160 
9 
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? 
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3 

? 
? 



(46) 
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(47) 



Question marks have been put into the entries which we do not know, but which do not 
affect the quantities calculated below. 

Transforming two-loop anomalous dimensions to the "standard" basis requires knowledge 
of one- loop diagrams with insertions of the extra evanescent operators Ek^\ ...,Eq'\ which 
introduces traces with 75 to the calculation. Fortunately, from among all the relevant quan- 
tities listed above, only cq^ and 034 depend on such traces. (In the cases of C54 and C74, the 
relevant color factors vanish.) The trace which occurs in the case of 054 is Tr(7^ 7^,75) which 
can be safely set to zero. On the other hand, 034 is irrelevant in the particular transformation 
we perform below. 
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At this point, we are ready to start the exphcit transformation. In the first step, we 
subtract some evanescent operators from the operators P5 and Pq 



(1) 
5 ) 



P' 



(1) 
6 ) 



(4J 



which is equivalent to taking the 6x8 matrix W in eqn. (|35| ) equal to 

k = 6, 



ik 



— 1 when i = k = 5 or i 
otherwise. 



(49) 



In the second step, we redefine the evanescent operators according to eqn. 



with 
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(50) 



Finally, we perform e- independent linear transformations according to eqn. (^if ), with 



and 



R 



M 



2 i 
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12 
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10 
10 
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0-^0 







-1 
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(51) 



(52) 



It is only a matter of tedious but simple algebra to verify that the triply primed basis is 



identical to the "standard" one, i.e. P!" = O, and E: 



,///(i) 



Of for i = 1, 6. The operators 
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Ej'^^^ and E^'^^^ play role of extra (unnecessary) evanescent operators in the triply primed 
basis, i.e. the physical operators do not mix into them at one loop. 

Finding the anomalous dimension matrix in the " standard" basis is now only a matter of 
simple matrix multiplication. Combining eqns. (0), and (^2p, one finds 



7' 



,///(0) 



7 



///(I) 



Af=ib + Wd 



where 



The formula (|54D resembles eqn. (3.4) of ref. 
of 7''^'* was considered in general. 

Substituting the matrices from eqns. (|2 
the above equations, one obtains 

-2 6 



Wc. 



(53) 
(54) 

(55) 



where renormalization-scheme dependence 
(H), (0) and dH) to 
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22 
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-2 + 1/ 


--f 


-f 

3'' 
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-6 


--f 


-f 


--f 


-16 + 1/ 



(56) 



and 
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Z+ 2 
2^3 
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65 


7 
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9 


3 




7 + 2 
2^3 
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21 


-f 


202 


1354 


1192 


904 
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243 


81 


243 


81 














5911 1 71 r 
486 ' 9 J 


5983 _|_ 1 f 
162 ~'~ 3'' 


2384 71 f 
243 9 J 


1808 1 r 
81 3'' 












379 _|_ 56 /■ 
18 ~'~ 243'' 


91 I 808 f 
6 ' SI J 


130 502 f 
9 243'' 


14 I 646 f 
3 81 














61 r 
9 


--f 

3 J 


Zl + 61 . 


-99 + ^/ 














682 J- 
243-' 


106 f 
81 


225 I 1676 r 
2 243 


1343 1 1348 r 
6 81 



(57) 



The above anomalous dimension matrices are in agreement with those found in refs. |10|, ^ 
in the so-called NDR scheme with fully anticommuting 75 (see eg. eqns. (4.5) and (4.8) in 
ref. §). 

Matching conditions in the "standard" basis can be obtained from eqns. (|Tl|)-([T^), ac- 
cording to 



1 



An' 



C{Mw) + 0{a 



(58) 
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with Af given in eqn. ([55|). The above equation yields 
11a 



1 



Stt 



S/5 



11a. 



w I 



247r 



1 



247r 



+ (59) 



+ (60) 



which agrees with refs. B [TO, |TTj| as well (see eg. eqns. (5.5) and (5.6) in ref 



5. Summary 

The I AF| = 1 nonleptonic effective hamiltonian has been considered in a renormalization 
scheme which allows to consistently use fully anticommuting 75 at any number of loops, but 
at the leading order in the Fermi coupling Gp. It is possible in the specific operator basis we 
have introduced. Such a basis is particularly suited for performing multiloop calculations. 

We have evaluated two-loop anomalous dimensions of the relevant operators and one-loop 
matching conditions for their Wilson coefficients in the new basis. Renormalization group 
evolution of the coefficients has been studied numerically in the case Ai? = —t\S = 1. 
In the end, the two-loop anomalous dimension matrix found by us was transformed to the 
previously used operator basis, and the existing NLO results have been confirmed. 
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Appendix 

Here, we give the nonphysical operators relevant in our discussion. In the process of renor- 
malizing off-shell two-loop amplitudes with insertions of operators Pi, ...Pq, there arise many 
EOM-vanishing operators, i.e. operators which vanish by the QCD equations of motion. 
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However, the specific structure of only one of them 

N, = -isLT-^^bL)D,G^f^'' + P4 (61) 
9 

is relevant in finding one- and two- loop mixing of the operators given in eqn. (^). 

Introduction of A^^i in section 2 was the reason why P4 (and, consequently, P3, P5 and Pg) 
entered the considerations. However, if P4 was not introduced this way, P3, ...,Pq would be 
generated anyway as counterterms to box diagrams (like the one in fig. |^) with insertions of 
the first term in A^^i. This term would be then treated as a physical operator. Introducing 
the nonphysical operator A^^i makes considering such box diagrams unnecessary, because the 



whole A^^i cannot have physical counterterms |18, 19, 21 




Figure 7: Sample box diagram with insertion of the first term in A^i (denoted by the black 
square) . 

Another kind of nonphysical operators which we need to list here are the so-called evanes- 



(62) 



cent operators i.e. operators which algebraically vanish in 4 dimensions [|1^, |2( 

The one-loop diagrams shown in fig. ^ with insertions of Pi,...,P6 require counterterms 
proportional to Pi, Pe themselves and to the following evanescent operators: 

Ei'^ = {sLl^,^^,^^,T''CL){cLl^'r'l^'T%L) - I6P1, 
E? = {sLl^.,l,,l,,CL){cLl^'l^'l^'hL) - I6P2, 

^3'^ = (^l7w7m2 7m3 7m47m5&l) i:,{qi^'Y'Y'l^'Y'q) - 2OP5 + 64P3, 

At the two-loop level, we encounter four more evanescent operators 

Ef^ = {sLj^,jf,,jf,,jf,,jf,,T''cL){cLY'Y'Y'Y'Y''T%^ - 20E[^^ - 256Pi, 
Ei'^ = {sLl^.,l^.,l^.,l^.,l,,TcL){cLY'Y'Y-'Y'Y'hL) - 204') - 256P2, 

= (sL7m7/.27/.37/.47M5 7M6 7M7&L) E,(g7'^^7'^^7^^7^*7^^7'^''7'^^g) - 336P5 + 1280P 



= (sl7w7^27^37/.47m5 7m6 7m7^"&l) E,(g7^^7'^^7^^7^*7'^^7'^«7'^^T"g) - SSGPg + I28OP4. 

(63) 

The latter evanescent operators E{ , ...,El are necessary in renormalizing the two- loop 
amplitudes from which one determines next-to-leading anomalous dimensions. They are 
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also present in eqn. pl| ) which determines the matrix c. However, a lot can be changed in 
their particular structure without affecting the two-loop anomalous dimension matrix of the 
physical operators. For instance, adding "e x (physical operator)" to E^"* affects d^^ and c 
in eqn. (|T^, but leaves 7*^^^ unchanged. This is in contrary to what happens when such a 
redefinition is applied to the one- loop evanescent operators E^"* (see eqns. (^8]) and (|39|) ). 

In section 4, the transformation of the two-loop anomalous dimension matrix of the phys- 
ical operators has been successfully performed without specifying what two-loop evanescent 
operators were used in the "standard" basis. This is another illustration of insensitivity of 
two- loop results to the specific structure of E^, . These operators become more important 
at the three-loop level. 
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